Entanglement of two coupled molecules 
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Entangled states of two coupled polar molecules controlled by laser pulses are studied theoretically. 
Schmidt decomposition is used to measure the degree of entanglement. By varying the pulse shape 
of the applied laser, transition from regular to chaotic behavior may occur. This quantum feature 
is due to the inference of different phases. Moreover, the entanglement is also found to be enhanced 
by increasing the strength of the laser pulses. 
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The control of the molecular alignment and orien- 
tation may be important in stereodynamics, trap- 
ping molecules, molecular focusing, surface catalysis, and 
nanoscale design. The alignment is due to the anisotropic 
polarizability induced by nonresonant pulses. The pen- 
dular states- hybrid of field-free molecular eigenstates 
IH El -can be created by turning on a picosecond 
laser pulse adiabatically. For polar molecules, a femtosec- 
ond laser pulse is able to generate a field-free orientation 
0- Q • The molecule remains oriented after turn off 
the pulse. Experimentally, several methods concerning 
the femtosecond photodissociation spectroscopy and the 
ion imaging have been employed to measure the rotation 
of molecules H U 03 El 

Recently, many researches have been focused to gen- 
erate entanglement in quantum-optic and atomic sys- 
tems 0, because entangled states are fundamental 
in quantum information processing |l2l . For molecu- 
lar systems, most works are focused on the implementa- 
tions of quantum computation . In general, the qubits 
are encode d by the vibrational (rotational) modes of the 
molecules 17]. In order to realize the quantum logic 
gates experimentally, the optimal control theory is usu- 
ally applied to design the suitable laser pulses. Recently, 
coupled-rotor- model attracts much interest also yjj be- 
cause some physical properties such as dielectric response 
may display peculiar behaviors in the presence of dipole- 
dipole interaction However, studies on the entangle- 
ment of coupled molecules under strong laser fields are 
still lack until now pfj| . In this rapid communication, a 
novel method to create entanglement between two cou- 
pled identical polar molecules separated in a distance of 
tens of nanometers is proposed. Both molecules are sup- 
posed to be irradiated simultaneously by the ultra-short 
laser pulses. It is found that the entanglement induced 
by the dipole interaction can be enhanced by controlling 
the inter-molecule distance or the field strength of the 
laser pulse. Moreover, the asymmetrical laser shape is 
also found to influence strongly the degree of the entan- 
glement. 

Consider now two identical polar molecules separated 
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in a distance of R. The molecule system is irradiated 
by a series of ultrashort half-cycle laser pulses j^J . The 
Hamiltonian of the system can be written as 
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Hf-rn, 
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where L'j and -|j (= B) are the angular momentum op- 
erator and rotational constant, respectively. Udi P is the 
dipole interaction between two molecules: 



Udi 
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(2) 



where fl 1 and \X 2 are the dipole moments. The dipole 
moments of two molecules are assumed, for simplicity, 
to be identical, i.e. = fi 2 = [i- The field-molecule 
coupling Hf- m can thus be expressed as 



IjlE (t) cos 9i cos (u>t) . 



(3) 



where 6q and 82 are angles between dipole moments and 
laser field. The laser profile is assumed in Gaussian 



shape, 



where Eq is the field 



i.e. E(t) = E e « 
strength, to is the center of peak, and a is the pulse du- 
ration. The time-dependent Schrodinger equation can be 
solved by expanding the wave function \1/ in terms of a 
series of field-free spherical harmonic functions Y\ m (0, <j>) 
as 



E 

iimii 2 m2 



(t)Y hmi X ) Y-; 2ro2 (0 2 , ^ 2 ). 



(4) 

where (6i,4>i) an d (02, ^2) are the coordinates of first and 
second molecule respectively. The time-dependent coef- 
ficients ci l7ni i 27n2 (t) corresponding to the quantum num- 
bers (lx, mi; ?2, TII2) and can be determined by solving the 
Schrodinger equations numerically. In the above equa- 
tion, the inter-molecule separation R is assumed to be 
fixed so that the total wavefunction has no spatial depen- 
dence. The variation of R might be inevitable due to the 
influence of laser fields or inter-molecule vibrations. How- 
ever, recent experiments exhibited that the spacial reso- 
lution in tens of nanometers for two individual molecules 
hindered on a surface is actually possible [2^, [2j| ■ The 



2 



free orientation model can be easily generalized to hin- 
dered ones, in principle, by replacing the spherical har- 
monic functions with hindered wavefunctions as shown 
in our previous work |2^|. Thus, the essential physics 
discussed here should be realistic. 

We now focus our attention on the entanglement 
generated in this system. The coupled molecules can 
be expressed as a pure bipartite system as = 
E; 1 m 1 ; 2 m 2 c hmihm 3 (t) \Y hmi ) \Y hm2 )- The reduced den- 
sity operator for the first molecule is defined as 



fcnoll =Tr mol2 |*) 



(5) 



Following the procedure of Schmidt decomposition, the 
bases of molecule 1 is rotated to make the reduced density 
matrix p mol 1 to be diagonal. The entangled state can be 
represented by a biorthogonal expression with positive 
real coefficients as 



I*) 



x lm) mo \ i l^m) mo i 2 1 



(6) 



where A; m is the eigenvalue corresponding to 
l y i™)moi 1 I^Wmol 2- Tne measure of entanglement 
for the coupled molecules can be parametrized by von 
Neumann entropy 



Entropy = - ^ A im log n A, 



! ill • 



(7) 
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The parameters for numerical calculations are based on 
Nal molecule whose dipole moment is 9.2 debyes and 
rotational constant is 0.12 cm^ 1 in the ground state. 
The duration and frequency of the half-cycle pulse are 
set to 279 fs and 30 cm -1 , respectively. The center of 
the peak is 1200 fs and the initial condition is set as 
coooo (t = 0) = 1. The ratio of the magnitude of posi- 
tive peak and that of negative peak of the pulse is 5 : 1 . 
Fig. 1 shows the time-dependent entropy after one pulse 
passes through this system. For inter-distance R = 5 
x 10~ 8 m, the entropy increases slowly from zero to al- 
most finite value of 0.4. For R = 1.5 x 10~ 8 m, on the 
contrary, the entropy grows rapidly with the increasing 
of time because the dipole force is stronger. Notes that 
the entropy only varies within a finite range at long time 
regime. This indicates that the systems reaches a dy- 
namic equilibrium state even though the dipole force is 
still present. In the insets, the orientations (cos(9i) and 
(cos #2), which can be evaluated immediately after the 
coefficients ci irrll i 2m2 (t) are determined, also reflect the 
effect of the dipole interaction. For R = 5 x 10 -8 m, 
the behavior of both molecules is quite close to that of 
a free rotor 6], in which the rotational period is close 
to nh/B. As the two molecules get close enough (the in- 
set of Fig. 1(b)), both molecules orient disorderly, and 
the periodic behavior disappears. This is because the 
dipole interaction is increased as the distance between 
the molecules is decreased, and the energy exchange be- 
tween two molecules becomes more frequently. The reg- 
ular orientation caused by the laser pulse is inhibited by 
the mutual interaction. 
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FIG. 1: Time evolution of the entropy after applying single 
laser pulse for (a) R = 5 X 10~ 8 m and (b) R = 1.5 X 10" 8 
m. The insets in (a) and (b) represent the orientations of 
two molecules at different distances. The field strength is set 
equal to 3 x 10 7 V/m. 



Fig. 2 illustrates the time evolution of the entropy with 
different ratios in magnitude of the positive and negative 
peak value of the laser pulse as R is set equal to 1.5 x 
10~ 8 m. For the case of ratio 9 : 1, an irregular- like 
behavior is still obtained, but its time-averaged value is 

0. 51, which is larger than the averaged one (0.43) in Fig. 

1. If the ratio is set equal to 1 : 1, the entanglement shows 
a nearly periodic behavior with small averaged entropy. 
This result is very similar to the case of no laser limit. 
This implies that the entanglement depends sensitively 
on the shape of the laser pulse. The positive and negative 
parts of the laser pulse seem to interfere with each other. 
Moreover, the dipole force only establishes a periodic-like 
entropy 

The distributions of the first ten contributive eigenval- 
ues (A p ) to the entropy at short and long time regimes are 
shown in the insets of Fig. 2. In the case of 9 : 1 ratio, 
the eigenvalue Ai dominates the contributions at short 
time regime(t = 50 ps). However, the dominant contri- 
butions are distributed to several values as t — 2000 ps 
regime. This means in long time limit the system is in 
some sort of dynamic equilibrium, and entropy saturates 
to certain value. On the contrary, Ai always dominates 
the contributions for either short or long time regime in 
the case of symmetrical laser shape as shown on the lower 
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FIG. 2: Time evolution of the entropy after applying single 
laser pulse for different ratios in magnitudes of the positive 
and negative peak value of the laser pulse. The graphs show 
the irregular (periodic) behavior for ratio 9 : 1 (1 : 1). The 
inset : the first ten contributive eigenvalues A p at short time 
(t = 50 ps) and long time (t = 2000 ps). The field strength 
and inter-molecule separation are E = 3 x 10 7 V/m and 
R = 1.5 x 10" 8 m. 



inset of Fig. 2. From statistical point of view, this some- 
how explains the suppressed and regular behaviors of the 
entanglement (entropy) . 

Fig. 3 shows the time evolutions of the populations 
of the eigenstates for different ratio of pulse shapes. For 
1 : 1 ratio, the energy level (0, 0; 0, 0) is mostly populated 
as shown in the lower panel of Fig. 3. Like the ground 
state, the populations of the higher levels (the inset of 
Fig. 3) also show the periodic behavior. The reason is at- 
tributed to the symmetrical shape of the laser so that the 
populations of the higher levels are almost contributed 
by the dipole interaction. In this case, the magnitudes 
of the higher level populations are rather small such that 
the periodic evolution of the entropy is obtained. On the 
other hand, for 9 : 1 ratio the populations of the higher 
states represent different degrees of irregularity as shown 
in the upper panel of Fig. 3. This is because a single laser 
ulse can generate high populations in the excited states 
Energy transfer by means of (mediated) dipole in- 
teraction generates the irregular evolutions of the higher 
excited states, which result in a randomly time-varying 
entropy. 

By adjusting the laser parameters, one can vary the 
degrees of the entanglement. Fig. 4 illustrates the time 
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FIG. 3: Populations of the states (Ji, mi; h, TOa) for different 
ratios. Upper panel : (1, 0; 0, 0) (dashed line), (1, 0; 1, 0) (solid 
line), (2,0; 1,0) (blue dotted line). Lower panel : (0, 0; 0, 0) 
(red line), (1, 0; 1, 0) (solid line), (1, 1; 1, 1) (dotted line). The 
inset in the lower panel is the enlarged figure showing the 
states (1,0; 1,0) (solid line), (1,1; 1,1) (dotted line), respec- 
tively. 



evolution of the entropy for fixed inter-molecule separa- 
tion R = 1.5 X 1CP 8 m and laser ratio 5 : 1 for a doubling 
of the laser field strength (Eq = 6 x 10 7 V/m). As can 
be seen, an irregular-like behavior of the entropy is ob- 
tained, and its averaged value is larger in comparison 
with Fig. 1(b). This can be understood well by study- 
ing the relationship between the dipolar interaction and 
the field strength. If the effect of laser field overwhelms 
the dipole interaction, the populations are distributed to 
a wider range. Since most of the populations are also 
distributed more averagely in this case, the entropy from 
Schmidt decomposition is certainly larger. Another way 
to control the degree of entanglement in this system is 
to change the positive and negative ratios of the laser 
pulse. Inset of Fig. 4 shows the time-averaged entropy 
with respect to different ratios. We find that the entropy 
is more enhanced as the ratio is larger. This means that 
the high asymmetric shape of laser pulse can generate a 
larger entropy at the same field strength. 

A few remarks about the differences between our model 
and previous works on generating entangled states should 
be addressed here. In our model, we consider the rota- 
tional excited states instead of internal electronic states 
of the molecules. Second, the laser frequency in our work 
is tuned as possible as far-away from resonance, while 
conventional creation of entanglement depends on the 
resonant driving pulses. This means our work provides 
a wider range to select the laser frequency to generate 
entanglement. As for the effect of decoherence, our en- 
tanglement is formed by the excited rotational states, in- 
stead of the vibrational states, therefore, the decoherence 
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is dominated by photon emission even if the molecules are 
attached to the surface of a solid. 



In conclusion, we have studied the entanglement of two 
coupled polar molecules irradiated by an ultra-short laser 
pulse. Schmidt decomposition is used to generated the 
time evolution of the entropy. The degree of the entangle- 
ment is analyzed by varying the inter-molecule distance, 
symmetry of the laser shape, and field strength. Its de- 
pendence on the controlled parameters may be useful for 
the quantum information processing. 



FIG. 4: Time evolution of the entropy for fixed field strength 
i?o = 6 x 10 7 V/m and ratio 5 : 1. The inset : Dependence of 
the time-averaged entropy on the pulse shape for fixed field 
strength Eo — 3 x 10 7 V/m and inter- molecule separation 
R = 1.5 x 10~ 8 m 
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